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Abstract 

This short essay discusses the application of Bogoliubov's theory of superflu- 
idity in the context of quantum phase transitions. 
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The importance of N. N. Bogoliubov's ground-braking paper On the Theory of 
Superfluidity [|l]] in the development of an understanding of superfluidity cannot be 
underestimated [^. More than 50 years after the publication of this seminal work, it 
continues to play a dominant role in contemporary condensed matter physics. It there- 
fore seems appropriate on the occasion of commemorating Bogoliubov's 90s birthday 
to submit a short essay discussing a modern application of his theory in the context of 
quantum phase transitions. Some of the material presented here is more extensively 
discussed in the review [|||. Other recent reviews can be found in Refs. [Q, ^, ^]. 

Bogoliubov's theory of superfluidity starts with the Lagrangian 

/: = 0*[zao-eHV) + Mo]0-Ao|0|*, (1) 

where the complex scalar field describes the atoms of mass m constituting the 
liquid, i9o is the total energy operator, while e(— iV) = — V^/2m is the kinetic energy 
operator, and /xq the chemical potential. The last term with a positive coupling constant, 
Ao > 0, represents a weak repulsive contact interaction. 

The theory features a global U(l) symmetry, under which the matter field acquires 
an extra phase factor 4){x) e*"</)(a:), with a the transformation parameter. Depend- 
ing on the ground state, which is determined by the minimum of the potential energy, 
the symmetry can be realized in two different ways. When fiQ < 0, the ground state is 
at = 0, and the system is in the symmetrical state. As the chemical potential tends to 
zero, the theory becomes critical, and when /io > 0, the global U(l) symmetry is spon- 
taneously broken by a nontrivial ground state, given by |0p = ^q/2Xq. This quantity 
physically denotes the number density ho of particles residing in the Bose-Einstein 
condensate. 

The spectrum of the single-particle excitations in this state is given by the cele- 
brated Bogoliubov form [|l]], 

E{k) = Ve'(k)+2A*oe(k), (2) 

whose most important signature is that at low momentum it takes the phonon form 
£'(k) ~ a/ fJ^o/m |k| predicted by Landau. The spectrum was shown by Beliaev [Q] 
to remain gapless when one-loop quantum corrections are included. And this was 
subsequently proven to hold to all orders in perturbation theory by Hugenholtz and 
Pines [Q], meaning that the Bogoliubov theory describes a gapless mode. This mode is 
nothing but the Goldstone mode accompanying the spontaneous symmetry breakdown 
of the global U(l) symmetry, and is the only degree of freedom present in this state. In 
other words, the Bogoliubov theory is a phase-only theory. At zero temperature and in 
the absence of impurities, the phase field is governed by the effective Lagrangian M 



eff 



2 m 



(3) 



where h is the average particle number density of the system at rest characterized by 
a constant phase field ip{x) — const, and c is the sound velocity, which to a first 
approximation equals c — y/JIofrn. The phase rigidity in the spatial directions, i.e.. 
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the coefficient of i(Viy9)^, is seen to be given by n/m, while that in the temporal 
direction by the compressibility k, because 



^ = n-'n. (4) 



Both these rigidities are response functions. Since the chemical potential jjL is repre- 
sented in the effective theory (|]) by | ]To[ ] 

= ~d(iLp{x), (5) 

a single differentiation of the effective Lagrangian with respect to ji yields the particle 
number density n{x) — n — {n/mc'^)[doLp + (V(^)^/2m] of the system slowly varying 
in space and time 

while a second differentiation yields the compressibility 



92 £, 



n^K, (7) 



as required. It also follows from Eqs. (||) and that n and are canonically conjugate 
variables [ 11 1. The form of the effective theory especially the combination d^^p + 
(V(/3)^/2m in square brackets is dictated by Galilei invariance 1 12 1. In cases where this 



symmetry is explicitly broken, as in the presence of impurities and at finite temperature, 
we expect changes in the relative weights of the coefficients (see below). 

Another, for the further development of the theory of superfluidity [||], momentous 
observation made by Bogolibov was the so-called depletion of the condensate. He 
showed that even at the absolute zero of temperature not all the particles reside in the 
ground state, but [|l]] 



no 3 \ TT 

where we replaced the coupling constant with the s-channel scattering length a = 
m\/2TT [ p3[ , p4| . (Recall that fiQ denotes the density of particles in the condensate.) 
Due to the interparticle repulsion, particles are removed from the condensate and put 
in states of finite momentum. In a strongly interacting system like superfluid ''He, the 
depletion is such that no more than about 8% of the particles condense in the zero- 
momentum state [ [Tst . 

Despite the depletion of the condensate, the phase rigidity in the spatial directions 
was found in Eq. ( 3|) to be given at the absolute zero of temperature and in the absence 
of impurities by the total average particle number density n/m. Since this coefficient 
denotes the superfluid particle number density (divided by m^), all the particles — 
not just those residing in the condensate — ^participate in the superfluid motion [p^. 
This changes at finite temperature and also when impurities are included: Galilei in- 
variance is broken then and ps no longer equals mn. On the other hand, the phase 
rigidity in the temporal direction as well as the first term in the effective Lagrangian 
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stay the same. This is because relation (|]) remains true. In general we thus have as 
effective theory [ p^ 

£cff = -ndo^ - + l^^^do^)^ + ■■■■ (9) 

Up to this point we have not specified the external parameter which must be var- 
ied to tune the chemical potential to its critical value where the system undergoes a 
phase transition. In the conventional application of the Bogoliubov theory, the control 
parameter is the temperature T. The critical temperature Tc can be determined within 
the theory by calculating the finite-temperature effective potential and identifying the 
temperature at which the minimum starts to shift away from the origin. At the one-loop 
level, one finds M: 



Tc = n 



^ n -2/3 1 / ,, X 2/3 9 Ci -\ 



where in obtaining this result a high-temperature expansion has been used. This is jus- 
tified because the leading term is of the order A^'^/"^, which is large for weak-coupling. 
Equation ( |lO| ) expresses the critical temperature in terms of the chemical potential. 
From the experimental point of view, however, it is more realistic to have the particle 
number density as independent variable. One then finds instead [P]: 



= CO [na^) ' , (11) 



where we again replaced A with the scattering length a, cq = — |C(^)/C(i) ~ 2.82, 

7 = ^, and Tq = (21: jm) is the critical temperature of a free Bose gas 

(A = 0). It follows that the critical temperature is increased by the weak repulsive 
interaction. This is qualitatively different from the strongly interacting ^He system. A 
free gas with ^He parameters at vapor pressure would have a critical temperature of 
about 3.1 K, whereas liquid ^He becomes superfluid at the lower temperature of 2.2 
K. A similar picture emerges from path-integral Monte Carlo simulations carried out 
by Griiter, Ceperley, and Laloe [p^. They found that at low densities, corresponding 
to small a, the critical temperature is increased by the repulsive interaction, while at 
higher densities it is decreased. In the weak-coupling limit, they found numerically 
the same exponent 7 = 0.34 ± 0.03 as in Eq. (pi]), while the value of cq was found 
to be an order of magnitude smaller: cq = 0.34 ± 0.06. As argued by these authors, 
a moderate repulsive interaction suppresses density fluctuations, resulting in a more 
homogeneous system. This facilitates the formation of large so-called exchange rings 
necessary to form a Bose-Einstein condensate. These exchange rings, as they appear 
in Feynman's theory of Bose-Einstein condensation [|9|, consist of bosons which are 
cyclicly permuted in imaginary time (see Ref. for a recent account). At higher 
densities, the exchange is obstructed because due to the strong repulsive interaction it 
is more difficult for the particles to move. This leads to a lower critical temperature. 

We now turn to the main subject of this essay, and consider the quantum critical be- 
havior of the Bogoliubov theory first studied by Uzunov [pT||. The critical behavior of a 



3 



system close to a quantum phase transition is dominated not by thermal fluctuations as 
in a classical phase transition at finite temperature, but by quantum fluctuations. In this 
context, the Bogoliubov theory is considered to be a phenomenological theory similar 
to the Landau theory of classical phase transitions. The system undergoes a quantum 
transition at the absolute zero of temperature when the chemical potential approaches a 
critical values /ic, which is not necessarily zero as in the case of the finite-temperature 
classical transition. The fine tuning of the chemical potential can be achieved by vary- 
ing a number of external parameters, such as the charge carrier density, the applied 
magnetic field, or the impurity strength. For values of the renormalized parameter 
larger than the critical value /i > Hc, the global U(l) symmetry is spontaneously bro- 
ken and the system is superfluid with a single-particle spectrum given by the gapless 
Bogoliubov spectrum, implying that the system is compressible. On lowering /i, this 
state is destroyed and replaced by an insulating state ^1%. 

In the absence of impurities, the insulating state is a so-called Mott-insulator, char- 
acterized by the absence of phase rigidity in both spatial and temporal directions, and 
by an energy gap in the single-particle spectrum. This insulating state, which arises 
solely due to the repulsive interaction, is consequently incompressible. 

On the other hand, in the presence of impurities, the bosons become trapped by the 
impurities, i.e., Anderson localized. The resulting insulating state is a so-called Rose 
glass characterized by a single-particle spectrum that is — like in the superfluid state — 
gapless. This state is therefore also compressible, so that the compressibility remains 
finite at the transition. 

To account for (quenched) impurities, the following term is added to the Bogoli- 
ubov theory: 

/:a = V(x)|(/'(x)p, (12) 
with '0(x) a real random field whose distribution is assumed to be Gaussian 



P{iIj) = exp 
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(13) 



and characterized by the impurity strength Aq. Physically, ip describes impurities ran- 
domly distributed in space. These impurities lead to an additional depletion of the 
condensate given in d space dimensions by [E3l E4] 



2<i/2-5^-<i/2p(2 _ d/2)m''/2A'^/2-2-d/2-i^^ 



The superfluid and normal mass density ps and respectively now become at the 
absolute zero of temperature | 



Ps = m\^n- ^"-aJ , Pn = -mnA- (15) 

It follows that the normal density is a factor A/d larger than the mass density mn^ 
knocked out of the condensate by the impurities. (For d = 3 this gives the factor | 
first found in Ref. [|2^.) As argued by Huang and Meng [p5|], this implies that part of 
the zero-momentum states belongs (for d < 4) not to the condensate, but to the normal 
fluid. Being trapped by the impurities, this fraction of the zero-momentum states are 
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localized. In other words, the phenomenon of Anderson localization can be accounted 
for in the Bogoliubov theory of superfluidity by including a random field. 

The universality class defined by the zero-temperature Bogoliubov theory is not 
only relevant to describe the critical behavior of superfluid films (either with or without 
impurities), but also to describe that of other systems, including Josephson junction ar- 
rays and superconducting films. In the so-called composite-boson limit, where Cooper 
pairs form tightly bound states, the BCS theory directly maps onto the Bogoliubov the- 
ory 1^], which is as we argued a phase-only theory. But even a weakly interacting 
BCS system was argued to be in the same universality class [|2^]. The reason is that the 
amplitude fluctuations of the order parameter are not critical at the transition, not even 
in the classical superconductor-to-normal transition in rf = 3 [^8|], only the phase fluc- 
tuations are. The phase of the order parameter therefore constitutes the relevant degree 
of freedom, which is precisely the one described by the Bogoliubov theory. (See, how- 
ever, Ref. [^], where it is argued that the amplitude fluctuations cannot be neglected, 
when considering quantum phase transitions in impure superconducting films.) The 
Bogoliubov theory presumably also forms the basis for the description of the critical 
behavior of fractional quantized Hall systems [^. 

To investigate the role of quantum fluctuations in the Bogoliubov theory we start 
with a dimensional analysis. Since, as far as the quantum critical behavior of this 
theory is concerned, the mass m is an irrelevant parameter, it can be scaled away by 
introducing t' = t/m, Hq = mji, Ag = Agm. The engineering dimension of the 
various variables is then easily determined as: 

[x] = -l, M = -2, M = 2, [Ao]=2-rf, [0] = id, (16) 

with d the number of space dimensions, and where we dropped the primes again. Note 
that the time dimension counts double as compared to the space dimensions. This is 
typical for nonrelativistic theories where the time derivative is accompanied by two 
space derivatives [see Eq. (|l])]. In two space dimensions, the coupling constant Aq has 
zero engineering dimension, showing that the |(/)|^-term is a marginal operator, and that 
dc = 2 is the upper critical space dimension above which the quantum critical behavior 
of the Bogoliubov theory becomes Gaussian. For d > quantum fluctuations are 
irrelevant, while for d < dc these fluctuations become crucial. 
Let us next compute the one-loop effective potential 

with i?(k) the gapless Bogoliubov spectrum (^. The integral over the loop momentum 
yields close to the upper critical dimension d — 2: 

Veff = -^--^^+0(e°), (18) 
4Ao ine k'^ 

where e — 2 — d, and k is an arbitrary renormalization group scale parameter, with the 
dimension of an inverse length. The rig ht-hand side of Eq. (|l|) is seen to diverge when 
the upper critical dimension is approached. The theory can be rendered ultraviolet 
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finite by introducing a renormalized coupling constant A 

1 K*^ m 

j = - + -, (19) 

where A = X/k'^. Its definition is such that for arbitrary d, A has the same engineering 
dimension as Aq in the upper critical dimension d — 2. As renormalization prescription 
we used the modified minimal subtraction. The beta function /3(A) follows as 



^eA+-A2. (20) 



Ao 



In the upper critical dimension, this yields only one fixed point, viz. the infrared-stable 
(IR) fixed point A* = 0. Below d = 2, this fixed point point is shifted to A* = en jm, 
implying that the system undergoes a 2nd-order quantum phase transition. Above the 
upper critical dimension, there is no (nontrivial) renormalization of the coupling con- 
stant, which explains why we omitted the subscript on /i and A in Eq. (|lo|). 



Since Eq. ( |18| ) could be rendered finite solely by a renormalization of the cou- 
pling constant, it follows that the chemical potential is not renormalized to this order. 
As shown by Uzunov these results remain true to all orders in perturbation theory [|T]]. 
The reason for this behavior is the special analytic structure of the nonrelativistic propa- 
gator at criticality, representing only particles propagating forward in time. As a result, 
the self-energy (and consequently ^) is not renormalized and the full 4-point vertex 
function is given by a geometric series, leading to the same beta function ( ^0| ) found at 
the one-loop order. Closely connected to this is that, despite the nontrivialness of the 
IR fixed point in c? < 2, the critical indices characterizing it are Gaussian [|2l|]. This 
conclusion was confirmed by numerical simulations in d = 1 [^1]]. 

This changes when impurities are included. A direct application of the renormal- 
ization group lead to the conclusion that the IR fixed point becomes instable. A 
more careful analysis, using a so-called double epsilon expansion, shows that the fixed 
point remains stable upon including impurities. The double epsilon expansion was 
originally introduced in statistical mechanics by Dorogovtsev [ ^3| ] to treat impurities 
of finite extend in a classical system. To consistently account for these in perturbation 
theory, one must assume their dimensionality to be small, and perform in addition 
to the usual epsilon expansion, also an expansion in e^. The impurities described by 
Eq. (|l2|) are static grains which trace out a straight worldlines when time is included. 
In other words, the impurities are line-like in spacetime, and have also to be treated in 
a double epsilon expansion, assuming that their dimensionality ed is not 1, but small 
instead. The quantum critical behavior of the Bogoliubov theory in d space dimensions 
with randomly distributed static impurities tracing out "worldlines" of dimensionality 
Cd falls in the universality class of a d-dimensional classical system with randomly dis- 
tributed extended impurities of dimensionality 2ed — at least to the one-loop order [^. 
The factor 2 arises because, as we mentioned before, in the nonrelativistic Bogoliubov 
theory, time dimensions count double as compared to space dimensions. 

Besides having a diverging correlation length ^, 2nd-order quantum phase transi- 
tions also have a diverging correlation time ^t, indicating the time period over which 
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the system fluctuates coherently. The way the diverging correlation time scales with 
the diverging correlation length, 

(21) 

defines the so-called dynamic exponent z. The traditional scaling theory of classical 
2nd-order phase transitions is easily extended so as to include the time dimension [p2|. 
Let (5 cx if — Kc, with K the external control parameter, denote the distance from 
the phase transition, so that ^ ^ ^ with v the correlation length exponent. At the 
absolute zero of temperature, a physical observable 0{kQ, |k|, K) at finite energy fco 
and momentum k can in the critical region be written as 

0(fco,|k|,if)=e'°0(6fco,e|k|), (r = 0), (22) 

where do is the scaling dimension of the observable O. The right-hand side does 
depend not explicitly on K, but only implicitly through ^ and 

Since a physical system is always at some finite temperature, we have to investigate 
how the scaling law ( ^2| ) changes when the temperature becomes nonzero. The easiest 
way to include temperature in a quantum field theory is to go over to imaginary time 
T = it, with r restricted to the interval < t < /3, where (3=1 /T is the inverse 
temperature. The time dimension thus becomes compactified. The critical behavior 
of a phase transition at finite temperature is still controlled by the quantum critical 
point provided < /3, so that the system does not notice the finite extend of the time 
dimension. Instead of the zero-temperature scaling (p^, we now have the finite-size 
scaling 

0{ko,\KK,f3)=P''°/'0{pko,P^/'\k\,p/^t), (T^O). (23) 

The distance to the quantum critical point is measured by the ratio ~ |(5|^''/r. 

Let us apply these general considerations to the effective theory (^ |l7|| . The 
singular part of the free energy density /sing, which scales near the transition as 

/sing~r^'+^\ (24) 

arises from the low-energy, long-wavelength fluctuations of the Goldstone field. The 
ensemble averages give 

ii^v?) - r', {{d,vf) - er' - r'"- (25) 

Combined, these hyperscaling arguments yield the following scaling of the rigidity 
constants: 

The first conclusion is consistent with the universal jump in the superfluid density pre- 
dicted by Nelson and Kosterlitz [|5|| for a KosterUtz-Thouless phase transition which 
corresponds to taking z = and d = 2. 

In an impure system undergoing an Anderson transition, the compressibility n'^n is 
nonsingular at the critical point and hence z ^ d for repulsively interacting bosons in 
an impure media [|4|. Surprisingly, the same conclusion holds for an impme fermionic 
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system [^j]. For d = 1 it follows that space and time appear symmetric as in a relativistic 
theory. 

In a clean system, on the other hand, with a density-driven Mott transition, i.e., 
5 oc /i — /ic, /sing can also be directly differentiated with respect to the chemical 
potential to yield for the singular part of the compressibility 

n^AJsing^ (27) 



In this case n'^K ^ n^Ksing, so that zv — \ |17| in accord with the Gaussian values 
= z = 2 found by Uzunov |21 1 for the pure case in d < 2. 



The above hyperscaling arguments have been extended by Fisher, Grinstein, and 
Girvin to include a l/|x|-Coulomb potential. This potential is important for the 
quantum phase transitions in charged systems because the Coulomb repulsion sup- 
presses fluctuations in the charge density and simultaneously enhances those in the 
canonically conjugate variable 0, thereby disordering the ordered state. The quadratic 
terms of the effective theory in Fourier space become when the l/|x|-Coulomb poten- 
tial is included 

4ff - \ (p^k^ - l^kl Ikl'^-i) |^(fco,k)p, (28) 

where e is the renormalized charge. Using similar hyperscaling arguments as before, 
one finds that this charge scales as 

- ^1-^ (29) 
Arguing that in the presence of random impurities the charge is nonsingular at the 



transition, the authors of Ref. [36| concluded that 

z = 1. (30) 

This again is an exact result which replaces the value z = dof the neutral system in an 
impure media. 

Most experiments on quantum phase transitions in charged systems measure the 
conductivity a. To describe such type of systems, we minimally couple the Bogoliubov 
theory to an electromagnetic vector potential [Aq, A). The conductivity turns out to be 
related to the superfluid mass density via | ^ ] 

(31) 

On account of the scaling relation (p^), it then follows that 

ar-^C'^'^'^K (32) 

implying that the conductivity and therefore the resistivity is a marginal operator in two 
space dimensions [0. 

The magnetic field H scales with f as iJ ~ $o/^^, where $0 = 27r/e is the 
magnetic flux quantum. This implies that the scaling dimension of A is unity, 

dA = 1, (33) 
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so that I A| ^ ^ ^. From this it in turn follows that the electric field £^ = |E| scales as 
E ^ Cr^'?^^ ^ £,^^^^^K and that the scaling dimension (Iao of Aq is z, 

d-Ao = z, (34) 

so that ^i"^ ^ 

Let us now be specific and consider quantum phase transitions triggered by chang- 
ing either the applied magnetic field, i.e., 6 cx iJ — i/c, or the charge carrier density, i.e., 
6 (x n~ Uc- For DC (A;o = 0) conductivities in the presence of an external electric field 
E we have on account of the general finite-size scaling form ( p3| ) with fcg = |k| = 0: 

a{K, T, E) = <?(|5r7r, W'^'+^^/E). (35) 

This shows that conductivity measurements close to a quantum critical point of the kind 
discussed here should in general collapse onto two branches when plotted as function 
of the dimensionless combinations |(5|'^^/T and |(5|''^^+^-' / E: a lower branch bending 
down for the insulating state and an upper branch tending to infinity for the other state. 
The best collapse of the data determines the values of vz and v{z + 1). In other words, 
the temperature and electric-field dependence determine the critical exponents v and z 
independently. 

The table below shows experimental data for the critical exponents z and v of the 
superconductor-to-insulator transition in thin films, the Hall-liquid-to-insulator transi- 
tion in fractional quantized Hall systems, and the conductor-to-insulator transition in 
silicon MOSFET's at extremely low electron number densities. 



Transition 


z 


V 


Superconductor-to-Insulator [ 




1.0 ±0.1 


1.36 ±0.05 


HaU-Liquid-to-Insulator |^ 


m 


w 1.0 


w 2.3 


Conductor-to-Insulator [p3 


n 


0.8 ±0.1 


1.5 ±0.1 



A few remarks are in order. First, the values for the dynamic exponent z found in 
these systems are in accordance with the prediction z = 1 recorded in Eq. (^0|), which 
was obtained using general hyperscaling arguments for an impure system with a l/|x|- 
Coulomb potential. Second, the values of the critical exponents characterizing the 
Hall-liquid-to-insulator transition are universal and independent of the filling factor — 
whether an integer or a fraction. Third, earlier experiments on siUcon MOSFET's at 
lower densities seemed to confirm the general believe, based on the work by Abrahams 
et al. [p4|], that such two-dimensional electron systems do not undergo a quantum 
phase transition. In that paper, where electron-electron interactions were ignored, it 
was demonstrated that impurities always localize the electrons at the absolute zero of 
temperature, thus excluding conducting behavior. Apparently, the situation changes 
drastically at low electron number densities, where the l/|x|-Coulomb interaction be- 
comes important. The values of the critical exponents found for this transition are 
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surprisingly close to those found for the superconductor-to-insulator transition. Since 
further experiments in an applied magnetic field p^ ] also revealed a behavior closely 
resembling that near the superconductor-to-insulator transition, it is speculated that the 
conducting state in silicon MOSFET's is in fact superconducting. 
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